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Abstract 

We test the disk/corona model of Janiuk et al. (2000) against the composite mean 
optical/UV/soft X-ray spectrum of radio loud and radio quiet quasars from Laor et 
al. (1997) which applies to faint quasars. The model well represents the optical/UV 
continuum if the hardening of the locally emitted spectrum is included, with the 
color to effective temperature ratio / ~ 2 in the inner 20-30 Schwarzschild radii. 
Comptonization seen in the soft X-ray band is well explained by the adopted corona 
model in the case of radio loud objects. However, this Comptonization is much 
stronger in radio quiet objects and additional ad hoc assumed Comptonizing medium 
must be present. We speculate that perhaps in both types of quasars there is a strong 
outflow of the hot plasma; this plasma is collimated in radio loud objects but not 
collimated in the radio quiet objects so it is present along the line of sight, serving 
as the required Comptonizing medium. Hard X-ray power law is not explained by 
the model and it may come from non-thermal component of electron plasma, as it 
is the case for galactic black holes in their soft states. 

Key words: Accretion, accretion disks. Galactic nuclei. Active galaxies. Quasars, 
Radiation mechanisms 
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1 Introduction 

The standard model of central en- 
gine of active galactic nuclei (AGN) 
consists of a massive black hole, an 
accretion disk, and a distant molecu- 
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lar torus. A major problem, however, 
is the location of the source of X-ray 
emission not expected from this sim- 
ple picture. X-ray emission carries a 
significant fraction of the bolometric 
luminosity. In Seyfert 1 galaxies it is 
mostly in the form of a hard X-ray 
power law component with the high 
energy cut-off at about 100 keV. In 
Narrow Line Seyfert 1 galaxies and 
quasars the hard X-ray power law 



Preprint submitted to New Astronomy 



1 February 2008 



contains much lower fraction of the 
bolometric luminosity but there is a 
second component in those sources - 
large soft X-ray excess - which car- 
ries large fraction of the luminosity. 
This X-ray emission comes from a 
relatively compact region located 
not too far from the accretion disk, 
as suggested by the fast variability 
in X-ray band and the presence of X- 
ray reflection component. However, 
even the extensive monitoring cam- 
paigns did not solve yet the problem 
of exact geometry of this region and, 
subsequently, of the physical process 
leading to the formation of the hot 
plasma. 

In the present paper we consider a 
particular model of the origin of this 
X-ray emission which is based on 
physical grounds. It is a model of a 
hot dissipative corona above an ac- 
cretion disk, as described by Janiuk 
& Czerny (2000). This model is based 
on assumptions that the corona is 
in hydrostatic and radiative equilib- 
rium with the disk, the rate of energy 
dissipation in the corona as well as in 
the disk is proportional to the pres- 
sure (total pressure in the disk and 
gas pressure in the corona) and the 
spontaneous division of the flow into 
disk and corona reflects the thermal 
instability within the irradiated gas, 
as discussed by [Krolik, IVlcKee ^ 



Tarter (1981|) . This model predicts 
the overall shape of the spectrum 
optical/UV/X-ray from an accretion 
flow as a function of the mass of the 
black hole, an accretion rate and the 
viscosity parameter a. We introduce 
some modifications to the model and 
we compare it with the observational 
data. 



One of the source of difficulties in 
modeling AGN is the fact that signif- 
icant fraction of the bolometric lumi- 
nosity is emitted in the unobserved 
region between far UV and soft X- 
rays so the bolometric luminosity 
and the spectral shape of the prin- 
cipal component - Big Blue Bump 
- are poorly determined for a sin- 
gle object. Therefore, in the present 
paper we decided to test the model 
against composite quasar spectra of 
Laor et al. (1997| ) which practically 
bridge the unobserved gap due to 
combining the spectra coming from 
objects at various redshifts. 

In Section 2 we shortly describe the 
composite and a theoretical model. 
Section 3 describes our results. We 
discuss the results in Section 4. The 
conclusions are given in Section 5. 



2 Method 



2. 1 Composite spectrum 



The mean spectra of radio quiet and 
radio loud quasars were shown by 
Laor et al. (1997 ). They cover a broad 
frequency range log v ~ 14.3 — 18.5, 
i.e. they extend from 15 000 A up to 
10 keV. 



The mean X-ray spectrum contains 
the spectra of objects from the Bright 
Quasar Survey ( phmidt fc Green 



with redshift z < 0.40 and col- 
umn density iV^f < 1.9 x 10^° cm~^. 
It was obtained adopting the fol- 
lowing cosmological parameters: 
the Hubble constant equal 



2 



5U kms ^IVlpc \ the deceleration Czerny fc Zycki, 1997|; |Janiuk k~Cz 



parameter q = 0.5, and the cosmo- erny, 20001 ). The total flux generated 



logical constant A equal to = 0.0. 
These low redshift objects deter- 
mine the mean spectrum above 



log V 



16.7 (200 eV). Below" 



it we have the unobserved range 
due to the Galactic absorption 
which extends from log ^ 15 to 
logv ^ 16.7. This range, however, is 
partially covered from 2000 A to 350 
A {log V ~ 15.2 — 16.0) by spectrum 
from [Zheng et al. (199"7| ) . This result 
is based on observations of quasars 
at high redshifts [z > 0.33). 

Final spectrum still contains a rela- 
tively short uncovered range. How- 
ever, the simple extrapolation of the 
spectra from both sides roughly coin- 
cide which strongly suggest that the 
resulting spectra are a good repre- 
sentation of the broad band spectra 
of bright radio quiet and radio loud 
AGN. 



2.2 Disk/ corona/ extended medium 
structure 



in the corona, F^, and in the disk, 
Fd, is determined by the standard 
formula (see e.g. [Kato, I'ukue fc 
neshige, 19981) 



Fd + F, 



3GMM 



g{r) 



where G{f) represents the New- 
tonian boundary condition at the 
marginally stable orbit 



Q{r) = 1 



/3Rs, 



chw 



(2) 



2.2.1 Local disk spectrum 

We assume that the disk is relatively 
cold, geometrically thin and optically 
thick (i. e. a-disk theory, [Shakura"^ 



Sunyaev, 1973| ) . We approximate the 
distribution of the angular momen- 
tum in the disk and the corona by 
Keplerian distribution 



IGM 



K 



(3) 



We consider a model of stationary 
accretion flow onto a Schwarzschild, 
massive black hole which consists 
of an accretion disk, a hot corona 
and (optionally) a spherically sym- 
metric extended hot medium. The 
accretion is ultimately responsible 
for the emission of the radiation, and 
the energy is dissipated partially in 
the disk and partially in the corona. 
Therefore, the accretion rate, M, to- 
gether with the mass of the central 
black hole, M, are the basic param- 
eters of our model (see also |Wit"t^ 



where G is the gravitational con- 
stant, M - the mass of the black hole 
and r is radial radius of the object. 
The disk is the source of the soft pho- 
tons, Faofti which are Comptonized 
locally in the corona and afterwards, 
at a greater distance, by the spheri- 
cal cloud. The flux Fgoft consists of a 
flux Fd generated in the disk at any 
given radius is denoted by Fd and of 
a fraction of the flux F^. generated in 
the corona and subsequently inter- 
cepted and thermalised in the disk. 
This flux determines the effective 



3 



temperature of the disk at a given 
radius 



soft 



(4) 



The spectral shape of this radi- 
ation is sometimes approximated 
clS cl blackbody. However, when an 
electron scattering dominate signifi- 
cantly over free-free and bound-free 
absorption this assumption is not 
justified. If the number of scattering 
events is large the saturated Comp- 
tonization effect is the most impor- 
tant and it shifts systematically the 
energy of photons while still preserv- 
ing a 'black body' spectral shape. 
Shimura fc 'lakahara (1995| ) showed 
that this is an important effect in 
AGN accretion disks. 

We introduce the spectral hardening 
factor / in such way that the flux 
can be described by the equation 
([Shimura fc 'I'akahara, 19951) 



(5) 



where T^ff is the effective tempera- 
ture, Tcoi = fTeff is the color tem- 
perature, and is the Planck func- 
tion of blackbody. 

The value of / should come in princi- 
ple from the radiative transfer com- 
putations. However, the results of 
Shimura & Takahara (1995) were 
obtained assuming pure hydrogen 
and helium atmosphere while heavy 
elements (particularly iron) are very 
important. Other codes 



le.g. 



Mer- 



loni, Fabian & Ross, 2000|;|Sobczak et 


al., 1999; 


Madej & Rozahska, 2000) 



mostly concerning the description 
of bound-free transitions and lines 
of heavy elements which are essen- 
tial. It is therefore more convenient 
to parameterize this effect, compare 
the results to the data and 'a pos- 
teriori' conclude which theoretical 
predictions are most appropriate. 

This effect is particularly important 
in the innermost part of the disk 
where the ionization is the highest 
but it should become negligible fur- 
ther out. Therefore we assume that 
the spectral hardening factor / is the 
monotically decreasing function. It 
decreases from the value fmax at the 
marginally stable orbit radius to the 
value / = 1 further away, at r larger 
than Vcoi- We use basically a linear 
function but results obtained with 
other functional shapes are similar 
(see Sect. 3). 



2.2.2 Local corona 

The corona is hot, optically thin 
and two-temperature (i. e. the ions 
temperature Tj is much greater than 
the electrons temperature T^, as dis- 
cussed by ^hapiro, Lightman fc Eard 



ley (1976| )). The two-temperature 
corona model was considered by 
a number of authors 



laar' 



IMarashi (1991 


); Witt, Czerny & 


IZycki (1997 


); 


Esin et al. (1997 


)). 



Here we follow the formulation of 
Janiuk fc Czerny (2000| ). 



also contain several approximations. 



We assume that the coronal plasma 
is an isothermal medium in the verti- 
cal coordinate, z, but non-isothermal 
toward the radial coordinate, r (i. e. 
g = and ^ ^ 0). The flux of the 
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is transported directly to this parti- 
cles ( INarayan et al., 1998 ). In this 
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case the flux Fc is proportioned to 
a- viscosity introduced by |Shakura"& 



Sunyaev (1973| ) and the total pres- 



sure P. 



VLxd I P(z) dz 



(6) 



Fig. 1. A schematic diagram of ra- 
diation fluxes distribution in the 
two-temperature corona model. F^i and 
Fc are the fluxes produced, respec- 
tively, in the disk and in the corona. 
Fsoft = Fii+r]Fc{l — a) is the flux of soft 
photons, which are Comptonized with- 



where zo is the vertical coordinate at 
the basis of the corona. In our calcu- 
late we neglect the coronal radiation 
pressure and the magnetic pressure, 
and we employ the results of the ver- 
tical integration obtained by [Witt J 



out disk, rj is a fraction of coronal pho- 
tons directed towards the disk, a is a 
disk albedo. 

soft photons, Fsoft, produced in the 
disk is partly Comptonized within 
the corona and it is partly trans- 
mitted throughout. The scattered 
Comptonized coronal flux, F^, is par- 
tially directed towards the disk, rjFc, 
and partially directed outwards to 
an observer, (1 — 'r])Fc, as a hard X- 
ray and gamma-ray radiation (see 
Fig. 1). X-ray photons irradiating 
the disk are partially absorbed and 
thermalised, contributing to Fsoft, 
and partiahy reflected by the disk 
surface, depending on the value of 
the X-ray albedo, a. We assume in 
our calculations that t] = 0.5 and 
a = 0.2. Such a formulation of the 
radiative disk/corona coupling was 
introduced by Paardt 



Czerny fc Zycki (1997| , Appendix D). 



At present 



^Qk^PoH^ 



(7) 



where Pq = P{zo) is the gas pressure 
at the basis of the corona, H is the 
pressure scale height of the corona 
given by the ions temperature, Tj, un- 
der the assumption of the hydrostatic 
equilibrium 



H 



1 



'kF 



Qk V 



(8) 



Such a formulation is independent 
from the physical mechanism of the 
corona heating since the scaling with 
pressure, H, may correspond either 
to the accretion heating by viscos- 



ity (IZycki et al., 1995| ; [Witt, Czerny 



TTMiy: 



iviarasmi & Zycki, 1997|) or to the magnetic 



I'he local dux i^c is generated by a 
heating mainly of the ions, because a 
gravitational energy by accreting gas 



heating (e.g. [Svensson fc Zdziarski, 
T99§. 



The energy of ions in the corona is 
immediately transported to electrons 
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through the electron-ion Coulomb in- 
teraction. The ions are cooled and 
this cooling is described by the equa- 
tion ([Shapiro, Lightman fc Eardley, 
T97|) 



The ionization stage of the medium 
is conveniently expressed through a 
ionization parameter S. It is defined 
ratio of the ionization radiation 
pressure to the gas pressure 



^ ^ 3 k{T, - Te 
2 rriH 



1 + 



UeipoH— , (9) 



cPn 



'111 



where i/ei = 2AAxlO^^ poT ~^-^\nA[s-^] 
is the electron-ion coupling rate and 
InA «i 20 is the Coulomb log. Here 
the density pn is related to the corona 



thickness and its optical depth Te. 



The heated electrons deliver the en- 
ergy subsequently to photons in the 
Inverse Compton process. In this way 
the part of observe/UV/soft X-ray 
photons emerge to a observator as 
a hard X-ray radiation. We describe 
this process as 



F, = A{Tes,Te,T,)K 



soft 



(10) 



where A{Tes,Te,Ts) is the Compton 
amplification factor and T, is the 
temperature of disk produced the lo- 
cal Fsoft. We interpolate the factor A 
from the tables prepared in advance 
from Monte Carlo computations. 

The vertical division of the medium 
into hot corona and relatively cold 
disk at every radius should not be 
arbitrary. Indeed, such a division re- 
sults naturally from the criterion of 
thermal instability in a irradiated 
medium studied by [Krolik, McKee fc 

whirh naturally HpvpI- 



I 'arfpr 



ops in the surface layers of an irra - 
diated accretion disk in hydrostatic 
equilibrium ( [Rozahska fc Czerny,| 



The phase transition from cold to 
hot medium occurs with a transcend 
of a specific value of the ionization 
parameter S ([Begelman, McKee &1 
Shields, 1983| ). 



E = 0.65 



10^ 



3/2 



(12) 



It scales with the electron tempera- 
ture and the fraction 0.65 is the value 
of the ionization parameter for an in- 
verse Compton temperature of 10^ K. 

The set of equations (3) - (12) al- 
lows to calculate the properties of the 
corona at each radius, function 
of global model parameters, i.e. mass 
of the central black hole, M, the ac- 
cretion rate, M, and the value of vis- 
cosity parameter, a, in the corona in- 
cluding the electron temperature and 
the optical depth. 

The local spectrum is now computed 
as a disk spectrum (Eq. 5) Comp- 
tonized by the local corona. Comp- 
tonization process was calculated us- 
ing code based on the semi-analytical 
formulae of [Czerny fc Z byszewska 



(1991| ) for a slab geometry. This de- 



scription allows for fast computa- 
tions on a fine radial grid. 
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2.2.3 Global disk/corona spectrum 



The final spectrum is calculated 
through the integration of the local 
spectra over all disk radii. Since the 
effect of corona is strongly dependent 
on the radius - it is weak close to 
the black hole and reaches maximum 
at a distance about a few hundred 
R-Schw, depending on the accretion 
rate, this integration is essential. It 
also means that hard X-ray emission 
comes mostly from the outer parts 
while weakest Comptonization in the 
innermost part creates a soft X-ray 
tail to the local disk spectrum. 



Extended 
medium 




2.2.4 Additional 
cloud 



Comptonizing 



Our corona model is universal in a 
sense that it does not contain any dif- 
ferentiation between radio loud and 
radio quiet objects while the compos- 
ite spectra for the two classes of ob- 
jects is clearly ditterent. 



Fig. 2. A schematic diagram of place 
the additional Comptonize cloud de- 
note as the extended medium. This 
cloud has a spherical geometry. It's de- 
scribed by the Thomson scattering TTh 
and the temperature of electrons, Tem, 
this cloud. 

In this case they are visible outward 
as the additional X-ray radiation. 

In this case the Comptonization pro- 
cess is computed once per model 
so we use an accurate Monte Carlo 
method. We applied the numer- 
ical program, which was written 
by pierlinski (20001) and which is 
based on procedures by [Corecki H 
;WiIczewsErj™]). 



To model such systematic difference 
we eventually allow for the existence 
of an additional extended medium 
(see Fig. 2) with parameters not de- 
termined by the accretion rate. This 
medium may be connected with the 
observed warm absorber, as a hot 
medium confining cooler absorb- 
ing clouds (lOgle et al., 2000|) . This 
cloud has a spherical geometry but 
we consider also a slab geometry. 
Physically it is characterized by the 
Thomson opacity, T^xt, and the elec- 
trons temperature. Text- The corona 
and mainly the disk are the source of 
soft photons, which are scaterred by 
the electrons in extended medium. 



Such an accurate method would be 
too time consuming to apply to the 
Comptonization by the corona at 
each radius and it was not really 
necessary since the final integration 
over the disk surface (an a range of 
temperatures and optical depths) 
smeared any specific Comptoniza- 
tion features even if present in the 
spectra (mostly first scattering ef- 
fect). 



2.2.5 Hard X-ray radiation 

It is generally believed that the hard 
X-ray radiation in radio-loud objects 
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comes from the jet. The extension of 
the hard X-ray power law into high 
energies is not well constrained in 
quasars so the nature of this emission 
is not known. The studies of galactic 
sources in their soft states, however, 
suggest that in such a spectral state 
of an accreting black hole the hard 
X-ray emission comes from the non- 
thermal tail of electron distribution 
(e.g. Pierlihski et al., 1999|, for Cyg 
X-1). 

Therefore, we do not expect that our 
thermal corona model will explain 
the hard X-ray radiation. We treat 
this hard X-ray power law emission 
as an additional component and we 
add it to the resulting spectrum 
adopting the slope and normaliza- 
tion appropriate for radio loud and 
radio quiet objects separately. 



3 Results 



3. 1 Without extended medium 



We first discuss the basic version 
of the model, i.e. disk plus corona, 
which is described by the mass of 
the black hole, M, accretion rate, 
M, coronal viscosity parameter, a, 
and two parameters determining the 
color temperature distribution in the 
disk: f.max and r^oi- 

The UV part of the composite spec- 
trum is mostly determined by the 
first two parameters which govern 
the total luminosity and the posi- 
tion of the peak on vF^ diagram. 
However, the last two parameters 




Fig. 3. The hardening factor / as a func- 
tion of the disk radius r. 

are also important since the mul- 
ticolor blackbody is too narrow to 
disk represent the |Laor et al. (1997|) 
spectrum adequately. 

However, if we allow for the presence 
of the saturated Comptonization skin 
in the upper layers of the disk we can 
represent the spectra adequately. 

We describe the radial dependence 
of the color temperature to the ef- 
fective temperature ratio through a 
hnear function of value fmax at the 
marginally stable orbit, dropping to 
1 at Tcoi (see Fig. 3). 

Both parameters, fmax and Tcoi 
broaden the spectral distribution in 
the far UV range but the deviation 
of the shape is somewhat different in 
those two cases. In Fig. 4 we show 
the dependence of the spectrum on 
the maximum value of the color to 
effective temperature ratio, fmaxi 
and on the extension of the disk skin, 
^coi (upper and lower part, corre- 
spondingly). We also experimented 
with other functional dependences of 
the parameter / on the disk radius 
(see Fig. 3) but when the parame- 
ters were adjusted to model the data 
possibly well the differences in the 
final spectral shapes obtained with 
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Fig. 4. The dependence of the spec- 
tral shape on the adopted value of 
Tf-oi radius (top) for f^ax = 2.5 and 
on the value of jmax (bottom) for 
r^oi = 25 Rschw Other parameters: 
mass M = 1.4 X 10^ Mq, accretion rate 
M = 1.5 M0yr~-^, viscosity parameter 
a = 0.08. To guide the eye, we also 
show a composite spectrum for radio 



loud and radio quiet objects from Laor 
et al. (1997] ) . 



different /(r) profiles were negligi- 
ble. Therefore, in further discussion 
we use always the linear formula. 

The value of the viscosity parameter 
a in the corona has no influence on 
the optical/UV part of the spectrum 
but it is important for the relative 
efficiency of the X-ray radiation. We 
choose it to reproduce the observed 
level of the soft X-ray emission in the 
Laor et al. (1997|) spectrum possibly 
well. 
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Fig. 5. The best representation 
of the composite spectrum of ra- 
dio quiet objects from [Laor et al. 
(1997] ) (short-dashed line) with the 



disk/corona model (continuous line). 
Model parameters: M = 2.4 x lO'^ Mq, 
M = 1.83 Meyr-i, a = 0.17, 
fmax ^ 2.3, real = 23 Rschw For com- 
parison, we show the disk blackbody 
spectrum (dashed-dotted line) and 
disk blackbody-l-corona spectrum (long 
dashed line). 

The best representation of the ob- 
served composite spectrum for radio 
quiet objects is shown in Fig. 5. The 
representation of the data is satisfac- 
tory in the optical/UV band in the 
X-ray band, where both the normal- 
ization and the slope are well repro- 
duced. The required value of the fmax 
parameter implies the color to tem- 
perature ratio equal 1.85 at 10 Rschw 

However, there is a large discrepancy 
between the predicted shape and 
the observed shape in the soft X-ray 
band. 

In our model the soft X-ray emission 
comes from the Comptonization of 
the disk emission in the innermost 
part of the flow. Clearly, the coronal 
strength in this region (i.e. the value 
of the Compton parameter deter- 
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Fig. 6. The best representation 
of the composite spectrum of ra- 
dio loud objects from |Laor et al. 
(19971 ) (short-dashed line) with the 



disk/corona model (continuous line). 
Model parameters: M = 2.9 x 10^ Mq, 
M = 1.68 M0yr-^ a = 0.2, 
fmax = 2.2, r^oi = 22 Rschw For com- 
parison, we show the disk blackbody 
spectrum (dashed-dotted line) and 
disk blackbody -|-corona spectrum (long 
dashed line) 

mined by the temperature and the 
optical depth) is too low to model 
well the observations. Since our basic 
model does not contain additional 
free parameters apart from the global 
parameters listed at the beginning 
of Sect. 3.1 we have no way to in- 
crease this effect artificially. It means 
that additional Comptonization is 
required, possibly by the extended 
hot medium. 

However, the mean quasar spectrum 
of radio loud objects in the soft X- 
ray band is well reproduced by our 
basic model (see Fig. 6). It means 
that in this case no extended medium 
is needed to explain the broad band 
spectrum. 



3.2 Extended medium in radio quiet 
objects 



The model consisting of the disk 
and corona emission as described in 
Sect. 2.2 did not account for the ob- 
served shape of the composite spec- 
trum of radio quiet quasars in the 
soft X-ray band. Now we consider a 
case when the disk/corona system is 
embedded in an extended hot opti- 
cally thin cloud. 

We fix the parameters of the disk / corona 
system at the previous values ap- 
propriate for radio quiet and radio 
loud objects, and we add two new 
parameters: the temperature of the 
extended medium, Te^t, and its opti- 
cal depth. Text- 



The soft X-ray emission in the com- 
posite spectrum of |Laor et al. (1997|) 
is not exactly of a power law shape 
but shows a discontinuity in the slope 
at logu ~ 16.7. Therefore, if we fol- 
low this shape precisely we can de- 
termine both parameters of the ex- 
tended medium independently. This 
would not be the case if the emission 
is a single power law since in such case 
only the slope, and subsequently the 
combination of the two parameters, is 
determined. It is difficult to say if this 
change of slope is real - it is caused 
by bridging the unobserved gap with 
a power law. However, the values of 
Text and Text stiU may be quite repre- 
sentative of the required medium. 



The result of this modeling for ra- 
dio quiet objects is shown in Fig. 7. 
Not surprisingly, this model well rep- 
resents the overall spectrum. The hot 
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Fig. 7. The best representation 
of the composite spectrum of ra- 
dio quiet objects from [Laor et al 



(19971 ) (short-dashed hne) with the 
disk/corona model plus extended 
medium (continuous line). Model 
parameters: M = 2.4 x 10^ Mq, 
M = 1.83 M0yr-\ a = 0.17, 

fmax — 2.3, T'col — 23 Rschwi 

Text = 60 - 75 keV, Text = 0.55 - 0.5. 

medium parameters: Text = 60 — 75 
keV and Text = 0.55 — 0.5, respec- 
tively. 

The exact values of the extended 
medium parameters depend on the 
adopted geometry. For a comparison, 
we also considered a flat configura- 
tion for the extended medium. The 
quality of the data representation is 
the same. The best values of the Text 
and Text are equal 55 keV and 0.45. 



4 Discussion 



We considered a physically justified 
model of the disk/corona. The emis- 
sion from such a system calculated 
in the simplest direct way did not 
account for the extended shape of 
the spectrum in far UV and for the 
shape of the soft X-ray spectrum in 



radio quiet and radio loud objects in 
the composite spectra of P^aor et al 
TT997I ). 



Two modifications to the simplest 
picture were necessary: (1) the pres- 
ence of the hotter disk skin where 
saturated Comptonization leads to 
a noticeable difference between the 
color temperature and the effective 
temperature at r < 30 Rschw (2) the 
presence of the hot extended medium 
surrounding the disk and the corona 
in the case of radio-quiet objects. 



4-1 The reliability of the composite 
spectra 



Our results depend sensitively on the 
spectral shape of the composite spec- 
tra adopted as the representation of 
a typical spectrum of radio loud and 
radio quiet quasar, respectively. 

The reliability of the composite spec- 
trum of [Laor et al. (1997| ) may be 
easily questioned since various types 
of objects were used to construct two 
major parts: high redshift quasars 
(z ~ 1) were used by Zheng et al. 
(1997) to obtain the far UV part 
and low redshift objects were used 
to find the soft X-ray part. This was 
unavoidable since this was the only 
way to close the observational gap 
between the two bands. Also the 
number of objects used to prepare 
the spectrum (~ 100) was not very 
large. 



The confidence in the composite may 
be further decreased by the presence 
of systematic differences between this 
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composite and other composites in 
the optical/UV band. 

In Fig. 8 we show the comparison 
of schematically drawn Zheng et al. 
(1997) composite and the optical 
part of the Laor et al. (1997) against 
the Francis et al. (1991; upper panel) 
and Morris (private communication; 
middle panel) . We see strong flatten- 
ing in Zheng et al. which lead to the 
conclusion about high color to effec- 
tive temperature ratio in quasar ac- 
cretion disks. Almost no such flatten- 
ing is present in two other compos- 
ites: standard black body accretion 
disk is almost a perfect representa- 
tion of the data for Francis et al. (see 
Fig. 17 in the review of Koratkar"^ 
Blaes, 19991 ). 



The difference, however, is most 
probably due to the difference in 
the studied population of quasars. 
Francis et al. (1991) sample con- 
sists of quasars at z ~ 3 and of the 
bolometric luminosity of order of 
10^^ — 10^^ erg s~^ and Morris com- 
posite is based on the same sample 
while |Zheng et al. (19971) part of the 
Laor composite was based on z ~ 1 
quasars of bolometric luminosity of 
order of 10'^^ erg s~^. 

As an argument in favor of this in- 
terpretation we plot in Fig. 8 the op- 
tical/UV spectrum of PG1211+143 - 
a radio quiet quasar classifled also as 
NLSl galaxy. The bolometric lumi- 
nosity of this object is 2 x 10^^ erg s~^ 



( Janiuk, (Jzerny fc Madejski, 20(IT ). 
The shape of this spectrum is notice- 
ably similar to the composite and 
cannot be modeled with a black body 
accretion disk. It is deflnitely flat- 




15 15.5 16 
log{u) [Hz] 

Fig. 8. The comparison of the Zheng et 
al. (1997) composite in the UV/far UV 
band (thick solid line) and Laor et al. 
(1997) in the optical band (thin solid 
line) with Francis et al. (1991) compos- 
ite (upper panel), Morris (private com- 
munication) composite (middle panel) 
and the spectrum of a quasar/NLSl 
PG1211 from Elvis et al. (1994, optical 
data) and Bechtold et al. (2001, HST 
data). 
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ter in UV than Francis et al. (1991 



composite. On the other hand, in- 



dividual spectra of high redshift, 
bright quasars show on average much 
harder spectra, with mean value of 
the slope of 0.36 (i.e. oc z/~o-36) for 
a total of 24 radio quiet quasars, con- 
sistent with [Francis et al. (1991| ), in 
the observed 1000 - 2000 A band in 
quasar rest frame ([Scott et al., 2000[ ) . 



The dependence of the spectral shape 
on the luminosity state is well known 
for galactic sources in their high/soft 
states dominated by the disk emis- 
sion. Sometimes the soft component 
is well described as disk black body 
emission but frequently additional 
Comptonization of this component is 
needed, independent from the hard 
X-ray tail (for a review see e.g. [Zycki,[ 
2001[ ). We therefore conclude that 
our results apply to faint or mod- 
erately bright quasars, and do not 
apply to bright quasars. 



Unfortunately, we cannot repeat the 
whole analysis for bright objects. The 
composite of [Francis et al. (1991|) 
appropriate for bright quasars is not 
supplemented by the X-ray part. 
The data on single bright objects 
suggest that the hard X-ray slope 
([Reeves fc 'I'urner, 2000|) is similar 



to that given by Laor et al. (photon 
index F = 1.89 ± 0.05 for radio quiet 
quasars and 1.66ib0.04 for radio loud 
quasars) . The relative level of hard X- 
ray emission in comparison to UV in 
very bright objects may also be sim- 
ilar: for example, very bright quasar 
PG 1247+268 has loguL^ equal to 
47.09 at 3000 A ( [INeugebauer et al. 



198^ ) and 45.72 at 2 keV ( [Reeves fc 



index as determined by [Laor et "aT 
(1997[ ). However, the soft X-ray part 



becomes unconstrained since the soft 
X-ray excess moved into unobserved 
range. Therefore, it cannot be used 
effectively to test the corona model. 



4- -2 Mean quasar properties 



Our disk/corona model implies that 
the studied quasar sample does not 
accrete at a rate close to the Edding- 
ton limit. The value of the mass of 
the representative black hole is about 
3 X 10^ Mq for radio loud objects and 
somewhat lower for radio quiet ob- 
jects (~ 2.5 X 10^ Mq), and the lumi- 
nosity to the Eddington luminosity 
ratio is equal 0.016 and 0.022, corre- 
spondingly, for a non-rotating black 
hole and pure hydrogen plasma. This 
ratio is rather low, typical for Seyfert 
1 galaxies, only the mass of the black 
hole and consequently, the luminos- 
ity is proportionally higher. This re- 
sult is a direct consequence of the 
high color to effective temperature 
ratio since roughly 



L/LEdd oc 



(13) 



'i'urner, 



meaning the same 



Disk blackbody fits would therefore 
indicate the ratios 3-4 times higher 
but we did not find them satisfactory 
in case of [Laor et al. (199"7[ ) com- 
posite. Our value of the black hole 
mass agrees with the value obtained 
by Kawaguchi et al. (2000) since the 
model contained similar color to ef- 
fective temperature ratio. Their ac- 
cretion rate for radio quiet quasars 
is slightly higher {L/LEdd = 0.031) 
because their model effectively ac- 
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counted also for the soft X-ray ex- 
cess, although at the expense of arbi- 
trary assumptions about the corona 
properties. 

Fits to Francis et al. (1991) compos- 
ite would give similar value of the 
black hole mass (see e.g. [KoratkaFS 
Blaes, 1999| ) but much higher accre- 
tion rate, and the L/LEdd ratio closer 
to 1. 

The value of the viscosity parameter 
determined from the model {a ~ 0.2) 
is reasonable. Similar values are de- 
rived for accretion disks in cata- 
clysmic variables (e.g. Smak, 1999) 
and in galactic X-ray novae during 
outburst (e.g. Lasota, 2001). Models 
of fast variability of the microquasar 
GRS 1915+105 under the influence of 
radiation pressure instability suggest 
somewhat smaller value, a ~ 0.01 
(Janiuk, Czerny & Siemiginowska, 
2000). Determinations available for 
AGN based on the observed optical 
variability supposed to correspond 
to thermal timescale give rather 
weak constraints [a > 0.01, Siemigi- 
nowska & Czerny, 1989; a ~ 0.03, 
Webb & Malkan, 2000). 

The value of the viscosity parameter 
is essential from the point of view 
of the inner boundary condition of 
an accretion disk. Low values, below 
0.1, are consistent with the standard 
description of the accretion in the 
vicinity of the marginally stable orbit 
while a > 0.1 leads to a considerable 
loss of angular momentum below the 
marginally stable orbit due to the 
type of the sonic point ( [Muchotrzeb 



field (Krolik, 1999; Reynolds, Ar- 
mitage & Chiang, 2001). Determined 
values are just at the border. 



4-3 Disk skin 



Large color to temperature ratio re- 
quired by the data is an interesting 
issue since it requires low level of ab- 
sorption and strong dominance of the 
disk atmosphere by electron scatter- 
ing. 

Significant ionization of the disk 
surface at small radii is not unex- 
pected. The observation of the X-ray 
reflected component is a complemen- 
tary way to estimate the disk opacity 
close to the surface observationally 
(for a review, see Mushotzky, Done & 
Pounds, 1993). Although the study 
of quasars are much more difficult 
than the study of Seyfert galaxies 
Nandra et al. (1997) and recently 
Reeves & Turner (2000) were able 
to determine the iron Ka line prop- 
erties on the basis of ASCA data. 
The line energy is at ~ 6.57 keV 
for objects of the brightness com- 
parable to those in the composite 
spectrum, which corresponds to the 
value of the ionization parameter ^ 
about a few hundred (e.g. Zycki & 
Czerny, 1994). The medium at the 
temperature ~ 10^ K and optical 
depth Tes ~ 1 is expected to develop 
under moderate external irradiation 
by the (non-thermal?) hard X-ray 
component since the temperature 
corresponds to an Inverse Compton 
temperature for effectively soft total 



Czerny, 1986| ; |Kato, Fukue fc Mi- 



neshige, 1998| ) or strong magnetic 



radiation field (e.g. |Czerny & Du- 


mont, 1998; 


Nayakshin, Kazanas & 
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Kallman, 2000). 



corona is in the hydrostatic equihb- 
rium with the disk (PR.6zahska et ah 



Jiowever, such a skin will result in 



19990 . 



very weak Comptonization effect, 
and in particular, it will not change 
the color temperature of the Big Blue 
Bump peak. To explain the color to 
effective temperature ratio of about 
2 we need a thick layer {r^s ~ 100) 
remaining at the color temperature 
due to some heating and due to cool- 
ing by saturate Comptonization. The 
absorption in this zone should be 
negligible - otherwise the tempera- 
ture would drop down to the effective 
temperature. 



The results of Hhimura & 'i'akahara 



( 1995|) as well as [Kawaguchi, Shimura 
fc Mineshige (20001) predict such a 
strong Comptonization because the 
computations do not contain bound- 
bound and bound-free transitions 
which are efficient coolants at the 
discussed temperatures, particularly 
oxygen, carbon and iron lines are im- 
portant. Computations of irradiated 
media including those transitions 
do not show the color to effective 
temperature ratio in the thermalised 
emission significantly different from 
1 (e.g. Madej & Rozahska, 2000; 
Nayakshin et al., 2000; Dumont et 
al., 2001). The required level of ion- 
ization for efficient Comptonization 
seems to be higher than naturally 
predicted in an irradiated medium. 
Similar problem was found by |Cz- 
erny fc Zycki (1994| ) when they tried 
to explain the universal position of 
the soft X-ray excess as the reflec- 
tion feature from the weakly ionized 
medium. The presence of the heavy 
corona further suppresses the ioniza- 
tion level at the disk surface if the 



It is possible that the effect is caused 
by significant generation of the en- 
ergy close to the disk surface. Such 
an enhanced generation was not 
included in the performed radia- 
tive transfer computations. It would 
mean that the vertical distribution 
of heat generation is not well param- 
eterized by the aP prescription and 
the magnetic field may be responsi- 
ble for such a global redistribution. 

On the other hand, the conclusion on 
the large color to effective tempera- 
ture ratio relies strongly on the shape 
of the composite spectrum in UV - 
far UV band which still has to be 
tested (see Sect. 4.1). 



4-4 Soft X-ray discrepancy for 
radio-quiet objects 



The adopted corona model ade- 
quately described the soft X-ray part 
of the composite spectrum for ra- 
dio loud objects. Low optical depth 
of the corona and a broad electron 
temperature range predicted by the 
model at different radii was adequate 
to give the observed effect of moder- 
ate Comptonization. 

However, in this model the Comp- 
tonization is definitely too weak to 
account for the soft X-ray spectrum 
of radio quiet objects. Similar mod- 
eling of the composite spectrum of 
radio quiet objects was performed 
by |Kawaguchi, Shimura fc Mineshige| 
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(2000) and in their case the modehng 
was successfuL This was caused by 
the fact that in their approach the 
fraction of the energy dissipated in 
the corona and the optical depth of 
the corona at a given radius were free 
parameters of their model. In our 
corona model the radial distribution 
of those quantities were calculated 
from the global parameters (black 
hole mass, accretion rate and viscos- 
ity parameter) since the spontaneous 
division of the flow into the disk and 
corona part results from the cooling 
properties of the material in the op- 
tically thin and optically thick part. 
As a result, our corona was never as 
thick as assumed by Kawaguchi et 
al., independently from the assumed 
viscosity parameter. 

Our corona model, adopted after Ja- 
niuk & Czerny (2000), contained sev- 
eral simplifications. However, recent 
progress in corona modeling does 
not seem to improve the situation. If 
the corona is assumed to be accret- 
ing, and the advection is included, 
the corona above a disk is even 
weaker than in our description ( |Ja- 



niuk, Zycki fc Czerny, 20001) . Better 



description of the disk/corona tran- 
sition and replacement of Eqs. 11- 



Although a theory did not say its last 
word on this subject as the forma- 
tion of the corona is a badly under- 
stood process it may suggest that ac- 
tually the disk corona is too weak 
to produce the soft X-ray spectrum 
of radio quiet quasars and the exis- 
tence of the additional Comptonizing 
medium may be necessary. 

The comparison of the radio loud 
objects and radio quiet objects may 
suggest that the outflow takes place 
in all types of AGN (for a review, see 
e.g. [Veilleux et al., 2000|) . In radio 
loud objects the outflow is coUimated 
in a form of jet, so the jet contributes 
to the hard X-ray emission but gen- 
erally it does not cover our line of 
sight to the most part of the accre- 
tion flow. The mechanism of this 
collimation and its distance from the 
black hole is still under debate (e.g. 
Sikora, 2001). In radio quiet objects 
the outflow perhaps is not collimated 
so the outflowing plasma is in our 
line of sight and Comptonizes the 
disk/corona emission. We show such 
a putative geometry in Fig. 9. 

Some kind of outflow is frequently 
predicted by the corona models or 
two-dimensional computations of 



12 with the disk evaporation con- the accretion flow (e.g. Begelman, 


dition due to the electron conduc- McKee & Shields, 1983; 


Vleyer & 


tion (Kozahska & Czerny, 2000a) Mcycr-Hotmcistcr, 1994 




Igumen- 


leads again to the prediction of still schev, 2000 


; Turolla & Dullemont, 


weaker corona for high accretion 2000 


). However, the estimate of the 



rates (Rozahska & Czerny, 2000b). 
On the other hand, the presence of 
the magnetic field may weaken the 
efficiency of the conduction and re- 
duce the condensation effect found 
in the innermost part of the disk 
accreting at a high rate. 



amount of outflowing mass and its 
physical conditions is difficult and no 
predictions differentiating between 
the radio loud and radio quiet ob- 
jects are available within the frame 
of this approach. The properties of 
the outflowing material may signif- 
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Sir- 



Kallman, 2000|) since the hot layer is 




\ \ I / / 



RQ 




Fig. 9. Possible character of an outflow 
in radio loud (RL) and radio quiet (RQ) 
quasars. 

icantly depend on the considered 
distance from the black hole as well 
as the inclination angle which fur- 
ther complicates the issue (e.g. Elvis, 
2000). 



optically thin due to thermal insta- 
bility in the heated/cooled plasma 
(Krolik et al., 1981; Rozahska & Cz- 
erny, 1996). 

High resolution observations were 
performed so far in the case of a 
few Seyfert 1 galaxies, and the re- 
sults were not definite: Branduardi- 
Raymond et al., (2000) claimed that 
broad oxygen, nitrogen and carbon 
emission lines, of similar shape to 
iron line, are present in the XMM- 
Newton spectra of MCG-6-30-15 and 
Mrk 766 while Lee et al. (2001) ar- 
gue that they see only the warm ab- 
sorber effect on the continuum in the 
Chandra spectrum of MCG-6-30-15. 
More detailed analysis of new XMM 
and Chandra data should resolve 
this issue. 



4-5 Observational tests - soft X-ray 
emission lines 



The dominant role of the Comp- 
tonization in formation of the soft 
X-ray spectra implies that no strong 
emission lines are expected in this 
band. On the other hand, models of 
coronal flares predict that a signifi- 
cant fraction of the disk is not cov- 
ered by the corona but it is strongly 
irradiated. Such irradiated internally 
not heated layers are the source of 
rich emission line spectrum in soft X- 
ray band (e.g. [Nayakshin, Kazanas & 



5 Conclusions 

• disk/corona model of Janiuk & 
Czerny (2000) based on thermal 
instability well describes the com- 
posite spectrum of radio loud faint 
quasars in the optical-soft X-ray 
band if the color to effective tem- 
perature ratio is close to 2 in the 
innermost part of the disk 

• adopted disk/corona model does 
not account for the soft X-ray spec- 
trum in radio quiet faint quasars 
and additional Comptonizing 
medium is required, possibly due 
to the uncoUimated outflowing hot 
plasma 
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